The role of occupied d states in the relaxation of hot electrons in Au 
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We present first-principles calculations of electron-electron scattering rates of low-energy electrons 
in Au. Our full band-structure calculations indicate that a major contribution from occupied d 
states participating in the screening of electron-electron interactions yields lifetimes of electrons in 
Au with energies of 1.0 — 3.0 eV above the Fermi level that are larger than those of electrons in a 
free-electron gas by a factor of ~ 4.5. This prediction is in agreement with a recent experimental 
study of ultrafast electron dynamics in Au(lll) films (J. Cao et al, Phys. Rev. B 58, 10948 
(1998)), where electron transport has been shown to play a minor role in the measured lifetimes of 
hot electrons in this material. 



Relaxation lifetimes of excited electrons in solids with 
energies below the vacuum level can be attributed to a va- 
riety of inelastic and elastic scattering mechanisms, such 
as electron-electron (e-e), electraiL-phonon (e-p), and 
electron-imperfection interactions.ErH Besides, when these 
so-called hot electrons are generated by absorption of an 
optical pulse, as occurs in the case of time-resolised two- 
photon photoemission (TR-2PPE) techniques ,Bu elec- 
tron transport provides an additional decay component 
to the photoexcited electron population. Since inelastic 
lifetimes of hot electrons become infinitely long as they 
approach the Fermi level, e-p scattering and the scat- 
tering by defects both play a key role in the relaxation 
process of electrons very near the Fermi level. However, 
in the case of hot electrons with energies larger than 
~ 0.5 — 1.0 eV above the Fermi level, e-e interactions 
yield inelastic lifetimes that are in the femtosecond time 
scale and they provide the main scattering mechanism. 

Experimental femtosecond time-resolved photoemis- 
sion studies of electron dynamics have been performed 
in a variety of solid surfaces ^To the role of e-e inelas- 
tic scattering and that of electron transport being diffi- 
cult to identify. However, recent TR-TPPE experiments 
in Au(lll) films with thicknesses ranging from 150 to 
3000 A have shown the relaxation from electron transport 
to be negligible and the hot-electron lifetime to be solely 
determined, at energies larger than ~ 0.5 — 1.0 eV abowa 
the Fermi level, by e-e inelastic scattering processes.tLa 
Hence, these measurements provide an excellent bench 
mark against which to investigate the importance of 
band-structure and many-body effects on electron dy- 
namics in solids. Also, ballistic electron emission spec- 
troscopy (BEES) has shown to be capable of determining 
hot-electron relaxation times in solid materials x3 

In this paper, we report first-principles calculations of 
the energy-dependent inelastic lifetime of hot electrons 



in Au. We follow the-jnany-body scheipe first developed 
by Quinn and FerrellO and by Ritchie,&3 but we now in- 
clude the full band structure of the solid. This approach 
has already been successfully incorporated in the descrip- 
tion of inelastic lifetimes of excited electrons in a vacieta 
of simple (Al, Mg, and Be) and noble (Cu) metals£ac2l 
A similar methodology, which is also-based on the GW 
approximation of many-body theory,EJ has been used by 
other, authors to evaluate hot-electron lifetimes in Al and 

cu.y 

Our many-body scheme for the calculation of inel, 
lifetimes in real solids has been described elsewhere] 
For periodic crystals, the basic equation for the decay 
rate r~ x of an electron in the state 4>Kni °f ener gy ^k,m is 



(we use atomic units throughout, i.e., 
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where Bj/(q + G) represent the coupling between the 
hot-electron initial and final states [see Eq. (5) of Ref. 
|l9| , and £q G/ (q, u>) are the Fourier coefficients of the in- 
verse dielectric function of the solid, whichj-we compute 
in the random-phase approximation (RPA) £3 G and G' 
represent reciprocal lattice vectors, the integration over 
q is extended over the first Brillouin zone (BZ), the sum 
over rif is extended over the band structure for each wave 
vector in the first BZ, to — E^ ni — i5k-q.ii/ is the energy 
transfer (0 < u> < E^ ni — Ep), and Ep is the Fermi 
energy. 

If one neglects crystalline local-field effects, Eq. (|l|) 
can be expressed as follows 
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r -l = J Lv-/ d |^/(q + G)| 2 Im[e G , G (g^)] 
* 2 %Jbz ^ |q + G| 2 kG, G (q,w)| 2 ' 
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Here, initial and final states of the hot electron enter 
through the coefficients -Bi/(q + G). The imaginary 
part of eG,G(q, w) represents a measure of the number of 
states available for the creation of an electron-hole pair 
involving a given momentum and energy transfer q + G 
and u), respectively, renormalized by the coupling be- 
tween electron and hole states. The factor |eG,G(q> w )| 2 
accounts for the dynamically screened e-e interaction. 

The decay rate of Eqs. ([!]) and (|J) depends on both the 
wave vector k and the band index ni of the initial Bloch 
state. Since measurements of hot-electron lifetimes are 
reported as a function of energy and the proper choice of 
the wave vector k and the band index rii is usually not 
easy to be determined, we define t^ 1 (E) as the average 
of r _1 (k, n) over all wave vectors and bands lying with 
the same energy in the irreducible wedge of the Brillouin 
zone (IBZ). 

If all one-electron Bloch states entering both the coeffi- 
cients Bif (q+G) and the dielectric function eG,c(q, w) 
were represented by plane waves and, at the same time, 
all energy bands were replaced by those of free elec- 
trons, then decay rates would not be direction-dependent 
and both Eqs. (|l|) and (||) would exactly coincide|-|With 
the GW scattering rate of hot electrons in a FEGn For 
hot electrons with energies very near the Fermi level 
(E ~ £ , i?)_this result leads, in the high-density limit 
(r s <-p-,l)jHj to the well-known formula of Quinn and 
FcrrcllliZl [see Eq. (7) of Ref. [n| which yields time in fs 

when E and Ep are expressed in eV]. The r s 5 ^ 2 scal- 
ing of hot-electron lifetimes predicted by this equation is 
found to be entirely originated in the density-dependent 
screening of e-e interactions [high electron densities yield 
a weaker interaction, ~ r s 3 ^ 2 ] and the density-dependent 
momentum of the hot electron [high electron densities 
yield a larger momentum and, therefore, a smaller num- 
ber of available transitions, ~ r~ ]. 

Gold is a noble metal with entirely filled 5d-like bands. 
In Fig. 1 we show the energy bands of this face-centered 
cubic crystal. We see large differences between the band 
structure of this material and that of free electrons. 
Slightly below the Fermi level, at E — Ep ~ 2 eV, we 
have d bands capable of holding 10 electrons per atom, 
the one remaining electron being in a free-electron-likc 
band below and above the d bands. These d bands, 
which are concentrated within a width of ~ 5.5 eV, are 
associated with the characteristic d-band wave functions 
at levels r' 25 and T^. Hence, a combined description 
of both delocalized 6s 1 and localized 5d 10 electrons is 
needed in order to address the actual electronic response 
of Au. The results presented below have been obtained 
by first expanding all one-electron Bloch states in a plane=. 
wave basis and then solving the Koiiti-Sham equationES 
of density-functional theory (DFT)Ea with all 6s 1 and 



5<i 10 Bloch states taken-as valence electrons in the pseu- 
dopotential generation.E3 Thus, a kinetic-energy cutoff as 
large as 75 Ry has been required, thereby keeping ~ 1200 
plane waves in the expansion of each Bloch state. Though 
all-electron schemes, such as the full-potential linearized 
augmented plane-wave (LAPW) method,Ej are expected 
to be better suited for the description of the response of 
localized d electrons, the plane- wave pseudopotential ap- 
proach has already been successfully incorporated in the 
description-af-thc dynamical responses and hot-electron 
lifetimesBBB of Cu. 

The sampling over the BZ required for the evaluation of 
both the dielectric matrix and the hot-electron decay rate 
of Eqs. ([!]) and (||) has been performed onal6xl6xl6 
Monkhorst-Pack (MP) mesh.E3 For hot-electron energies 
under study (E — Ep ~ 1.0 — 3.0 eV), well-converged re- 
sults have been found with the inclusion in the evaluation 
of the dielectric matrix of conduction bands up to a max- 
imum energy of 25 eV above the Fermi level. The sums 
in Eqs. (|l]) and (||) have been extended over 15 G vec- 
tors, the magnitude of the maximum momentum transfer 
q+G being well over the upper limit of ~ 2 q F (q F is 
the Fermi momentum). 

Our full band-structure calculation of the average life- 
time t{E) of hot electrons in Au, as obtained from Eq. 
(|l|) with full inclusion of crystalline local-field effects, is 
presented in Fig. 2 by solid circles. The lifetime of hot 
electrons in a FEG with the electron density equal to that 
of 6s 1 electrons in Au (r s = 3.01) is exhibited in the same 
figure, as obtained from either Eq. ([!]) or Eq. (||) with 
the full RPA dielectric function of free electrons (solid 
line) and from the formula of Quinn and Ferrell (dotted 
line). Also plotted in this figure is the lifetime of hot 
electrons in Au(lll) films, as determined from accurate 
TR-TPPE experiments (open squares) li3 The agreement 
between our full band-structure calculation and the ex- 
perimental data is excellent, for all hot-electron energies 
under study. 

Both band-structure (solid circles) and FEG (solid 
line) calculations presented in Fig. 2 have been carried 
out within the same many-body framework: The electron 
self-energy has been obtained in the GW approximation 
and on-the- energy- shell [deviations of the actual excita- 
tion energy from the independent-particle approximation 
have been neglected] , and the RPA has been used to com- 
pute the dielectric function of the solid. Hence, the ratio 
(~ 4.5) between our calculated ab initio lifetimes and 
the corresponding FEG calculations unambiguously es- 
tablishes the impact of the hand structure of the crystal 
on the hot-electron lifetimcEH Furthermore, since many- 
body effects beyond our GW-RPA scheme are expected 
to be within 20% of our calculations,!^ we conclude that 
large deviations of the experimental measurements from 
the FEG prediction are mainly due to band-structure ef- 
fects. Nevertheless, one must be cautious with the com- 
parison between our full band-structure calculations and 
the experiment. Since Au has a wide energy gap along 
the TL direction just over the Fermi level, the k\\ = pho- 
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toemission in Au(lll) cannot be associated, in the energy 
range under study (E — Ep ~ 1.0 — 3.0 eV), with hot- 
electron Bloch states in the same direction, and there- 
fore the proper choice of the wave vector k appears to be 
unclear. We have evaluated hot-electron lifetimes along 
various directions of the wave vector, and have found 
that differences between these results and the average 
lifetimes presented in Fig. 2 are within ~ 20%, which 
gives support to our comparison with the experiment £3 

In order to investigate the role that occupied d states 
plays in the relaxation of hot electrons in Au, now 
we neglect crystalline local-field effects and compute 
hot-electron lifetimes from Eq. (g). Scaled lifetimes 
t(E) x (E — Ep) 2 of hot electrons in Au, as obtained from 
Eq. (||) by replacing the hot-electron Bloch states and 
bands entering |-Bi/(q + G)| 2 by those of free electrons 
and the dielectric function in |eG.G(q, w )|~ 2 by that of a 
FEG with r s = 3.01, are represented in the inset of Fig. 
2 by open triangles. These calculated lifetimes, which 
have been obtained with full inclusion of the band struc- 
ture of the crystal in the evaluation of Im[eG,G(q 5 UJ )], are 
close to those obtained within the FEG model of the solid 
(solid line) , showing that the combined effect of the den- 
sity of states (DOS) available for e-h pair creation and the 
small overlap between d states below and sp states above 
the Fermi level nearly compensate. Hence, deviations 
of actual hot-electron lifetimes from FEG predictions are 
mainly originated in the deviation of the hot-electron mo- 
mentum from the FEG prediction [the momentum of hot 
electrons depends on the actual DOS at the Fermi level] 
and the participation of d electrons in the screening of 
e-e interactions. 

The effect of virtual intcrband transitions giving rise to 
additional screening is to largely increase ihe hot-electron 
lifetime in Au, as occurs in the case of CuJLj a noble metal 
with analogous electronic structure. Since the change 
8e\ in ei(w) due to the presence of d electrons in the 
noble metals_is_found to. be practically constant at low 
frequencies, Ej'E3 Quinr£3 treated the FEG as if it were 
embedded in a medium of dielectric constant eo = 1 + Se\ 
instead of unity. The corrected lifetime is then found 

1/2 

to be larger by roughly a factor of e , i.e., ~ 2.5 for 
both Cu and Au. Nevertheless, our first-principles evalu- 
ation of the wave- vector and frequency dependent dielec- 
tric matrix leads us to the conclusion that the role that 
occupied d states play in the screening of e-e interactions 
is much more important in Au than in Cu. This is an 
expected result, since Au 5d electrons lie further away 
from the nucleus than Cu 3d electrons, and 5d bands in 
Au are, therefore, more free-electron-like than 3d bands 
in Cu [dLbands in Cu are concentrated within a width of 
~ 3 eVjEZI smaller than in the case of Au] . 

The result we obtain with full inclusion of the band 
structure of the crystal in the evaluation of Eq. (g), 
but still neglecting crystalline local- field effects, is rep- 
resented in the inset of Fig. 2 by open circles. Our full 
band-structure calculation of Eq. (Q) for Au is repre- 



sented in the same figure by solid circles, showing that 
neglecting crystalline local-field corrections results in an 
overestimation of hot-electron lifetimes of 40% — 50%. 
These crystalline local-field effects, which partially com- 
pensate the pronounced participation of d electrons in the 
screened e-e interaction, are due to large electron-density 
variations in this material. 

In conclusion, we have presented full band-structure 
calculations of the inelastic lifetime of hot electrons in 
Au. We have found that a major contribution from oc- 
cupied d states participating in the screening of e-e in- 
teractions yields lifetimes of electrons with energies of 
1.0 — 3.0 eV above the Fermi level that are larger than 
those of electrons in a FEG by a factor of ~ 4.5, in agree- 
ment with the experiment. The effect of virtual interband 
transitions giving rise to additional screening does not de- 
pend on whether the hot electron can excite d electrons 
[the d-band scattering channel opens at ~ 2 eV below the 
Fermi level] or not. On the other hand, the actual density 
of occupied states available for real transitions is found 
not to play an important role in the relaxation mecha- 
nism. Consequently, actual lifetimes approximately scale 
as (E — Ep)~ 2 , as in the case of a FEG. That lifetimes 
of hot electrons in Au can be fitted within this scaling 
behaviour has also been demonstrated by a theoretical 
analysis of electron-electron mean free paths in BEES 
experiments.113 
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FIG. 2. Hot-electron lifetimes in Au. Solid circles represent 
our full ab initio calculation of t(E), as obtained after averag- 
ing r(k, n)^ 1 of Eq. (|l|) over wave vectors and over the band 
structure for each k. The solid line represents the lifetime of 
hot electrons in a FEG with r s = 3.01, as obtained within 
the full GW-RPA. Open squares represent the experimental 
measurements of Ref. [if]. The dotted line represents the pre- 
diction of Quinn and Ferrell. Differences between the FEG 
calculations represented by solid and dotted lines are due to 
the fact that the formula of Quinn and Ferrell is obtained for 
hot-electron energies very near the Fermi level (E ~ Ef) and 
in the high-density limit (r s — > 0). The inset exhibits scaled 
lifetimes of hot electrons in Au. Open circles represent our 
calculation of t(E), as obtained after averaging r(k,n) _1 of 
either Eq. (|]J) or Eq. (Q) over wave vectors and over the band 
structure for each k and neglecting crystalline local-field cor- 
rections. Open triangles represent the result obtained from 
Eq. (I2J) by replacing hot-electron Bloch states in Bij (q + G) 



and |e G ,G(q, c 



by those of free electrons, but with full in- 



clusion of the band structure of the crystal in the evaluation 
of Im[e G ,G(q, ui)]. 



FIG. 1. Calculated band structure of Au along certain sym- 
metry directions. 
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